Renormalization of the Spin-dependent WIMP scattering off nuclei 
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We study the amplitude for the spin-dependent WIMP scattering off nuclei by including the 
leading long-range two-body currents in the most important isovector contribution. We show that 
such effects are essentially independent of the target nucleus and, as a result, they can be treated as 
a mere renormalization of the effective nucleon cross section or, equivalently, of the corresponding 
effective coupling with values around 25%. 
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OO 1 !• INTRODUCTION 

The combined MAXIMA-1 Q, BOOMERANG 0, DASI @ and COBE/DMR Cosmic Microwave Background 
O ^ (CMB) observations [H imply that the Universe is flat [j| and that most of the matter in the Universe is Dark 
JL ■ i.e. exotic. These results have been confirmed and improved by the recent WMAP data Q- Combining the data one 
^"'l finds: 

i 1' n b = 0.0456 ±0.0015, n CDM = 0.228 ±0.013, n A = 0.726 ± 0.015 . 
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Since any "invisible" non exotic component cannot possibly exceed 40% of the above SIcdm j exotic (non baryonic) 
matter is required, i.e. there is room cold dark matter candidates or WIMPs (Weakly Interacting Massive Particles). 
Even though there exists firm indirect evidence for a halo of dark matter in galaxies from the observed rotational 
^-j- ' curves, see e.g. the review it is essential to directly detect such matter. The possibility of such detection, however, 
t-H . depends on the nature of the dark matter constituents and their interactions. 

Since the WIMP's are expected to be extremely non relativistic, with average kinetic energy (T) w 
50 keV(m WIMP /100 GeV), they are not likely to excite the nucleus, even if they are quite massive m WIM p > 100 
GeV. Therefore they can be directly detected mainly via the recoiling of a nucleus (A,Z) in elastic scattering. The 
t-H ' event rate for such a process can be computed from the following ingredients: i) An effective Lagrangian at the ele- 
mentary particle (quark) level obtained in the framework of the prevailing particle theory. The most popular scenario 
• *-H ' is in the context of supersymmetry. Here the dark matter candidate is the LSP (Lightest Supersymmetric Particle) 
In this case the effective Lagrangian is constructed as described, e.g., in Refs. 

[IM1- At least till supersym- 



metry is discovered in the laboratory (LHC), other approaches are also possible, such as, e.g., technibaryon [19l. |2C 
mirror matter [2lL [22[ , Kaluza- Klein models with universal extra dimensions [23l |24| ii) A well defined procedure for 
transforming the amplitude obtained using the previous effective Lagrangian, from the quar k to the nucleon level. 
To achieve this one needs a quark model for the nucleon, for the coherent mode see, e.g., 

HUMS!!]- Similar effects 

occur in the case of the spin induced matrix elements, iii) knowledge of the relevant nuclear matrix elements [28l [29|. 
obtained with as reliable as possible many body nuclear wave functions, iv) knowledge of the WIMP density in our 
vicinity and its velocity distribution. 

In the standard nuclear recoil experiments, first proposed more than 30 years ago [30(, one has to face the problem 
that the reaction of interest does not have a characteristic feature to distinguish it from the background. So for 
the expected low counting rates the background is a formidable problem. Some special features of the WIMP- 
nuclear interaction can be exploited to reduce the background problems, such as the modulation effect [3l|. [32H40j 
and backward-forward asymmetry ex pec ted in directional experiments, i.e. experiments in which the direction of the 
recoiling nucleus is also observed [4l| |42j - [53l | . 

There exists a plethora of direct dark matter ex per iments with the task of detec ting WIMP event rates for a va riety 
of targets such as those employed in XENON10 [EI^ENONIOOO^XMASS JHl , ZEPLIN (57|, PANDA-X [Hf, 
LUX j5<|, CDMS |6(|, CoGENT [Hj], EDELWEISS [62[, DAMA (gaUl, KIMS VM and PICASSO [H H3- 

For the interpretation of the experimental data, knowledge of the nuclear matrix elements is essential. In connection 
with nuclear structure aspects, in a series of calculations, e.g. in [38|, [68|,|69| and references therein, it has been shown 
that for the coherent contribution, due to the scalar interaction, the inclusion of the nuclear form factor is important, 
especially in the case of relatively heavy targets. They also showed that the nuclear spin cross sections, which may 
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be relevant for odd mass targets, are characterized by a single, i.e. essentially isospin independent, structure function 
and two static spin values, one for the proton and one for the neutron, which depend on the target. 

In the present paper we will focus on the evaluation of the spin induced matrix elements. In particular we will 
explore the currents for spin-dependent WIMP scattering off nuclei, in particular the most isovector channel. We 
will examine i) the modification of the rates at the one-body level by including a suitable nucleon form factor, which 
modifies somewhat the universality of the suitably normalized spin isospin structure functions and ii) include the 
leading long-range two-body currents, which are predicted in chiral EFT |70j. We will show that the two body 
effects lead to a quenching of the isospin one spin matrix elements. Since, however the resulting quenching factor is 
independent of the target, it can be absorbed into the corresponding nucleon cross section. We have not examined 
such effects on the isoscalar spin matrix element since this mode, if present at all, is suppressed due to the quark 
structure of the nucleon. 



II. THE DIFFERENTIAL EVENT RATE 



The expressions involving the differential are fairly well known (see e.g. [71| for a brief review in our notation). For 
the reader's convenience we will list the main features here. The event rate can be cast in the form: 
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where the first term represents the time averaged (non modulated) differential event rate, while the second gives the 
time dependent (modulated) one due to the motion of the Earth (see below). Furthermore for the coherent mode we 
get: 
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while for the spin dependent process we get: 
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In the above expressions /z r (/x p ) is the WIMP-nucleus (nucleon) reduced mass, A is the nuclear mass number, m x is 
the WIMP mass, p(x) is the WIMP density in our vicinity, assumed to be 0.3 GcV cm~ 3 , and mt the mass of the 
target. Furthermore one can show that 



„coh /|2 coh 

a A - A a N , 
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(4) 



with cr^ h the elementary nucleon cross section for the coherent mode. In the case of the spin the expressions for static 
nuclear cross section cr^" 11 gi yen below and: 



dt\ 
du J 



dh\ 
du J 



y|a 2 ^ii(w)*i(a v / ^), 



(5) 



The factor a/2/3 is nothing but vol \l (v 2 ) since in Eq. ([2]) J (v 2 ) appears. In the above expressions a — (y/^Urbvo)^ 1 , 
vo the velocity of the sun around the center of the galaxy and b the nuclear harmonic oscillator size parameter 
characterizing the nuclear wave function, u is the energy transfer E R in dimensionless units given by 



E R 
Qo(A) 



, Qo(A) = [mpAb 2 ]- 1 = 40 A^ 3 MeV 



(6) 



In the above expressions F(u) is the nuclear form factor and Fu(u) the isospin 1 spin response function. Note that 
the parameter a depends both on the WIMP mass, the target and the velocity distribution. Note also that for a given 
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energy transfer Er the quantity u depends on A. 

The functions ^o(x) and ^i(x) arise from the WIMP velocity distribution, via function g(v m i n , Ve(o)), which depends 
on the minimum WIMP velocity for a given energy transfer, i.e. 



/ A m p Er ,_, 

For the M-B distribution in the local frame it is defined as follows: 

g(v m in,v E (a)) = 1 a P"" e-^ +2V - v ^+^ a »K vdvd^ v max = v esc (8) 

{V^Vo) Jv min 

Ve(cx) is the velocity of the Earth, including the velocity of the Sun around the galaxy, Ve{ol) — eo(v s + 
S(smax — cosacos7y + cosa sin 71)5)) with 7 w 7r/6, £> s a unit vector in the Sun's direction of motion, x a unit 
vector radially out of the galaxy in our position and y = v s x x. S — 0.235 is the velocity of the earth around the sun 
divided vq and a is the phase of the Earth (a = 0, around June 3nd)*. The above upper cut off value in the M-B is 
usually put in by hand. Such a cut off comes in naturally, however, in the case of velocity distributions obtained from 



the halo WIMP mass density in the Eddington approach [72j, which, in certain models, resemble a M-B distribution 

To obtain Eqs ([4]) and ([5]) we expand g(v m i n , Ve{&)) in powers of 5, keeping terms up to linear in 8 w 0.135. We 
found it convenient to express all velocities in units of the Sun's velocity vq to obtain: 

vog(v m in,VE(ot)) = ^o(x) + *i(a;) cosa, x = (9) 

^o(x) represents the quantity relevant for the average rate and ^i(x), which is proportional to d, represents the effects 
of modulation. 

In the case of a M-B distribution these functions take the following form: 

ijf (x) = - [erf(l -x) + erf(x + 1) + erfc(l - y eso ) + eric{y CBC + 1) - 2] (10) 
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-erf(l — x) — crf(x + 1) — erfc(l — y csc ) — erfc(y eBC + 1) 

2 



'7T WTT WTT yJ-K 



(11) 



where erf(x) and erfc(x) are the error function and its complement, respectively, and y esc — v esc /vo, 550 < v esc < 660 
km/s. 



III. THE SPIN DEPENDENT WIMP-NUCLEUS SCATTERING 



The spin dependent WIMP-nucleus elastic scattering cross section can be derived from an interaction which at the 
quark level can be cast in the form 

£ = -A^J.^A g * ?775 * 9 (12) 
Where the scale parameter A as well as J depend on the particle model. Thus: 



1 One could, of course, make the time dependence of the rates due to the motion of the Earth more explicit by writing a fa 
(6/5)tt (2(t/T) - 1), where t/T is the fraction of the year. 
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1. Supersymmetry, see, e.g., a review in our notation [74| . 

A = 2(|C* 3 i| 2 -|C4i| 2 ), J = X775X (13) 

where C31, C41 are the Higgsino components of the WIMP, which in this case is the Lightest Supersymmeric 
Particle (LSP). The coefficients A q ca be read from the space component of the neutral current: 



J\ = -<ll\{\sin 2 e w - ~(1 - 75) - sin 2 9 w r 3 }q, A = 0, A x = ±, g° A (q) = 0, g\{q) = -1. (14) 
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2. Klauza Klein theories in models with Universal Extra Dimensions |24j |. 

In this case we encounter a number of cases depending on the nature of the WIMP, e,g. 

• The WIMP is a gauge boson. 
In this case: 

A « - tan 2 W — J = (*e x e , A u = —,A d = A s = —, A = -,A 1 = - 

3 \ m s / m s 18 18 2 3 

- = 1 (15) 

with m p the proton mass and e and e the polarization vectors of the initial and final WIMPs (the nor- 
malization factors 1 / y/rriB entering each boson field in the expression of the cross section has already been 
absorbed into A). 

• The WIMP is heavy neutrino. 

It is assumed the process proceeds via Z exchange. A viable possibility is if it is a Majorana Fermion. In 
this case: 

A = 1, J = —^775^ (16) 
The coefficients A q are given again by the space component of the neutral current (Aq = 0, A\ = |, g° A {q) = 

0, g\{q) = -1). 

In the present work we will not concern ourselves with any specific particle model, but we will, instead, focus on the 
following points: 

1. Going from the quark to the nucleon level, i.e. constructing an effective transition operator at the nucleon level. 
This involves: 

• The modification of the currents due to the structure of the nucleon. 

The effective couplings g\ (g A ) are obtained by multiplying the corresponding elementary amplitudes 
obtained at the quark level, as given above, by suitable renormalization factors g A and g A given in terms 
of the quantities Aq given by Ellis [zl,H, namely Am = 0.78 ± 0.02, Ad = -0.48 ± 0.02 and As = 
-0.15 ±0.02., i.e. 

9a(N) - g° A (q)(Au + Ad + As) - 0.13g° A (q), g\(N) = g\(q)(Au - Ad) = IMg^q) (17) 

The quantities g A (q) and g A (q) can be read from the hadronic current. For the neutral current they are 
unity, so the renormalization of the isovector component of the axial current is the usual one, while the 
isoscalar one is greatly suppressed, consistent with the EMC effect, i.e. the fact a tiny fraction of the spin 
of the quarks is coming from the spin of the quarks. It is for this reason that we started our discussion in 
the isospin basis. In the case of the boson WIMP in Kaluza- Klein theories one finds [24j |: 

g A (N) = ^Au+^(Ad + As) = 0.26, g A (N) = ^Au - ^Ad = 0.41 (18) 

Thus in general, barring very unusual circumstances at the quark level, the isovector component is expected 
to be dominant. 
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The nuclcon form factor 129, 76], which modifies the isovector axial current: 



9A q 
If we choose as a z-axis the direction of momentum we find: 



9A{q) ^ S = <T -i^l4 (i 9) 

_i_ /i-v-i Z 



1 q 2 +ml J TO — ^Q) 

[ cr m , m = ±1 

Thus only the longitudinal component of the transition operator is modified by the inclusion of the nucleon 
form factor. Thus, even at sufficiently high momentum transfers, only 1/3 of the differential rate will be 
affected. 

• possible effects arising from possible exchange currents, which lead to effective 2-body contributions. We 
will exploit recent work on Chiral Effective Field Theory (EFT) pTOj ] . 

2. Construct the relative Spin Structure functions. 

For this purpose and for illustration purposes we will consider the simplest possible nuclear system, which in 
our opinion is the A = 19 system [29| currently employed as a target by the PICASSO collaboration. The study 
of other experimentally more interesting targets , see e.g. ref. [66] ■ can be found elsewhere (70j . 

The nuclear ME entering the WIMP-nucleus cross section takes the form: 

\ME\ 2 = alS n (u) + ai a Soi(u) + a^S 00 (u), a = A g° A (N), a l = A x g 1 A {N), (21) 

i.e. the parameters ciq and a\ depend on the parameters A q as well as on the quark model for the nucleon. The 
structure functions 



with the functions Sy defined in the appendix, 2 are essentially independent of the isospin. Thus one can simplify the 
expression involving the nuclear ME entering the WIMP-nucleus cross section and write it as follows: 

\ME\ 2 = (affi? + aia O Oi + a 2 ,^ 2 ,) F n (u), (23) 

where f2o arL d fl\ are the static nuclear matrix elements involving cr for isospin and 1 respectively defined by 



n p = 2^ < \jj\s p \jj n n = 2^ t±l ^ \jj\s n \jj y, r» - a p + n n n 1 = n p - n n (24) 

As it is common practice the nuclear physics parameters f2 p and ft n (or fio and Oi ) as well Fu appear in the expression 
of the rate, while the parameters A, a\ and ao can be absorbed into the nucleon cross sections a p and a n (or ao(N) 
and cti(A)). Leaving aside the factor Fu, taken care of explicitly in Eq. ([5]), the static nuclear cross section becomes 



n\a x (N) + 2sign(a 1 a )n n 1 y/a (N)(n(N) + n 2 a (N) 



(25) 



with 



n 2 m 2 n 2 m 2 

<7o (N) = A 2 al^L, ai{N)= A 2 al^P^ (26) 

and mjv the nucleon mass. 

The structure functions Fij obtained in this work are shown in Fig. SJa), while those derived from a previous work 



2 The functions Sij for the Xe isotopes 1771 . [7811 are normalized differently. The functions Fij extracted from them are independent of this 
normalization. 
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|70| are exhibited in Fig. SJb) and (c). The isospin behavior of the functions JFy in the A = 129 is what we expect. 
We do not understand the stronger isospin dependence in the case of the A = 131. 

Before proceeding further we will interpret the results of the recent work (70j . summarized below (see section ITVT) . 
in our notation. If we define: 

r «<°> = §g (27) 

then 

1. for the 129 Xe isotope we take 

Q p = 0.034, fi„ = 1.140, n = 1.174, fti = -1.105, 

and 

roo(0) = r xl (0) = r O i(0) = ±- ~ 0.038 (1-body) , 

»7T 

roo(0) = ^-, r u (0)= 0.0240, 
^oi(0) = 0.0614 (one and two-body) 

2. for the 131 Xe isotope we have 

fl p = -0.023, 0„ = -0.702, Q Q = -0.776, = 0.679, 

and 

r 00 (0) = ru(0) = roi(0) = ±- ~ 0.079 (1-body) , 

47T 

roo(0) = ru(0) = 0.0487, r i(Q) = 0.1155 (one and two-body) 

47T 

The above results have been obtained from Eq. (16) as well as the data of Table 1 of Ref. [70j |. 

3. In the case of 19 F we have p9l |: 

VL P = 1.644, n n = -0.031, ft = 1-614, fii = 1.675 (28) 

The precise value of the parameters depends on the parameters of the model. For a reasonable choice of 
these parameters we get: 

roo(0) = r u (0) = roi(0) = 1 (1-body) , 

r 00 (0) = 1, ru(0) = 0.600, r i(Q) = 0.775 (one and two-body) (29) 

In the presence of effective two body terms we can write: 

rn(0)| (1+2)6 = (1 - <5) 2 r u (0)| 16 , r 01 (0)| (1+2)6 = (1 - 5) r 01 (0)| 16 (30) 

From the above analysis we find: 

5 « 0.229 (A=129), <5 w 0.234 (A=131), 5 w 0.235 (A=19) (31) 

The corrections due to the presence of two body terms are significant. From the above expressions, since 5 is 
independent of A, one has two possibilities: 
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1. Write Eq. in the form 



(oi (1 - 5fn\ + sign(a ai)yoT\/^)(l - 6)n Qi + cr ^o 



(32) 



where the nucleon parameters remain the same but the isovector nuclear spin matrix elements are quenched by 
a factor of (1 — S). 

2. View the correction factor like an effective charge. 

In this case, due to the nuclear medium, the elementary nucleon amplitudes are modified: 

ao — > a>Oi Oi— >oi(l — 5), (To ^ > °o, <ri — > eri(l — 5) 2 (33) 

but the nuclear spin matrix elements remain unchanged. In this case we recover the standard formula 

spi „ _ (gjjji + sign(aoai) v /aT v /g "»oOi + cr ^g) , 
°\4 — g 134J 

The parameter (5 depends on the assumed model as shown in table Q] 

One may express the above results in the proton neutron representation. The above behavior of the form factors 
imply that the proton and neutron form factors are almost the same. If one retains the amplitudes unchanged one 
finds: 
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(35) 



(36) 



In other words the two-body terms induce a small neutron (proton) component in those cases in which such a 
component is negligible in the standard treatment with only one body terms. 

On the other hand one can absorb the quenching into the coupling constants, in which case the proton and the neutron 
cross sections get modified : 



= a + ai(l - S), 



oq - ai(l - 5), cr p -> cr p (a + ai(l - <5)) 2 , o- n ->• <T„(a - ai(l - <5)) 5 



while the standard expression for the nuclear cross section is retained, i.e. 

°T n = n t^P^p + Vn^l + 2sign(a p a n ) y /opyfcr^ri p ri n ] 



(37) 
(38) 



IV. THE SPIN STRUCTURE FUNCTIONS IN THE CONTEXT OF EFT 



In this section we will employ the normalization of the structure functions even though, as shown above, they will 
be used only in extracting the relevant static spin matrix elements and in particular the quenching factors. Following 
the recent work of Ref. [7(| we have included the currents for spin-dependent WIMP scattering off nuclei at the 
one-body level as well as the leading long-range two-body currents, which are predicted in chiral EFT. The nuclear 
matrix elements (ME) entering the WIMP-nucleus cross section can be written as [76| 



\MEr = 



< j\m 



el 5 | 



J > 



< J 



\cl\ 



J > 



L = odd 



(39) 



where T£ 15 is the transverse electric and L\ the longitudinal projections of the axial-currents. Since the ground state 
of 19 F is the J = 1/2+ then only the L = 1 multipole contributes. As consequence the above operators can be written 
in the isospin basis in terms of the operators T^ L ' = ji{qri)\Y l {f\) x cr] L as 



1 1 

752 



1 



(40) 
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TABLE I: The parameter 8 entering the suppression factor / = (1 — ft) for the axial isovector contribution for various values of 
the model parameters p, I(p,P = 0), C3 and C4. 



p 


/ 


C3 


c 4 


5 


0.10 


0.58 


-3.2 


5.40 


0.256 


0.10 


0.58 


-2.2 


4.40 


0.203 


0.10 


0.58 


-3.40 


3.40 


0.189 


0.10 


0.58 
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2.40 
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0.10 


0.58 
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3.96 
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0.10 


0.58 


-3.78 


2.96 
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0.11 


0.59 
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5.40 


0.287 


0.11 


0.59 


-2.20 


4.40 


0.227 


0.11 


0.59 


-3.40 


3.40 


0.212 


0.11 


0.59 


-2.40 


2.40 
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0.11 


0.59 


-4.78 


3.96 


0.261 


0.11 


0.59 


-3.78 


2.96 


0.202 


0.12 


0.60 


-3.20 


5.40 


0.318 


0.12 


0.60 


-2.20 


4.40 
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0.12 


0.60 


-3.40 


3.40 


0.235 


0.12 


0.60 


-2.4 


2.40 


0.169 


0.12 


0.60 


-4.78 


3.96 


0.290 


0.12 


0.60 


-3.78 


2.96 


0.224 



and 



T£L\(q) = -j=l-J==J2 [ao(- T^) + V2T^) + an? h(q)(- + V2T^)] (41) 



where 



and 



ZgnpnK q 2 n p q 2 



2m N g A (mi + q 2 ) F£ 4m^ + q 1 



q 2 



/ 2 (g) = l-J-2-^. (43) 

The above operators are given in terms of the variable u = q 2 b 2 /2, where q is the momentum transfer and b the 
harmonic-oscillator length. Using the results of our previous shell model calculation for 19 F [2i| the parameter b has 
been taken the value 1.63 fm. The multi-particle matrix elements of the above operators in the considered sd model 
space are given in the Appendix. 

The resulting 2b contribution to the axial- vector WIMP current has been included as a density-dependent renor- 
malization 70] ai(l + 5), with 

^^^^ = 0)Q(2c 4 -c 3 ) + ^-). (44) 

A typical range of the densities in nuclei p — 0.10...0.12/m~ 3 has been taken into account. This leads to I(p, P) = 
0.64. ..0.66, using the Fermi-gas mean-value P 2 = 6k 2 /5 where kf the Fermi momentum, and I(p,P) — 0.58. ..0.60 
for P = 0. The coupling constants C3 and C4 describing the contributions of long-range one-pion-exchange as well as 
the short-range parts in the two-body currents are taken from Ref. In Fig. Q]the contour plot is shown for the 

parameter S as function of p and 2c 4 — C3. The curves in the plot correspond to constant values of 5. 

A second contribution <5af in the pseudo-scalar component of 2b current has been included via a density and 
momentum dependent modification: 

2 

fof (q 2 ) = 2c 3 -4 -r4 (45) 

U ' F 2 4m 2 . + q 2 v ; 
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FIG. 1: (Color online) Contour plot of the parameter S as a function of density p and the difference 2c4 — C3 of the parameters 
C3 and C4. The curves in the plot correspond to constant values of S and lighter shading denotes the increase of 5. 



where = 138.04 MeV and F v = 92.4 MeV. All the other parameters in Eqs. (|40l) and (|41l) of the appendix have 
been taken the values = 1040 MeV, g^pn = 13.05, qa = 1.27. The first ^-dependent term of /2(<?) and fi(q) arise 
arise at the one-body (lb) currents, while the second g 2 -dependent term in /2(<?) and the momentum independent 
term 5 are taken from the two body (2b) currents. All the terms of the form factors fi(q) and /2(g) are presented in 
Fig. [2] As it is seen the first <? 2 -dependent term of /2(g) is negligible. In Fig. [2] the corresponding average value of 
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the two body terms 6 and 6af is presented. At finite momentum p the structure functions Soo, Sii and Sqi can be 



TABLE II: Fits to the structure factors Soo, Su and 5*01 for spin-dependent WIMP elastic scattering off 1 F nucleus, including 
lb and 2b currents. The fitting functions are given in terms of an exponential e~" multiplied by a forth-order polynomial. The 
rows give the coefficients of the u" terms in the polynomial. 



rap 
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FIG. 2: (Color online) Form factors terms inserting in Eqs. (|42p and (|43p . 

given in terms of a forth-order polynomial multiplied by an exponential e~ u (see Appendix). The results of the fitting 
procedure are given in Table [TT1 while the structure functions Sij for lb and lb+2b currents are presented in Fig. [2] 
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FIG. 3: (Color online) Structure functions Soo,Sii and Soi for 19 F as a function of u — q 2 b 2 /2, for one body(lb) and one+two 
body(lb+2b) currents. No two body corrections were considered for Soo- 
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FIG. 4: (Color online) Structure factors Fqo, Fh, and Foi for 19 F, 129 Xe and 131 Xe isotopes and as a function of u = q 2 b 2 /2, 
for lb+2b currents. No two body corrections were considered for Fqo- 



DISCUSSION 



In the present paper we studied the evaluation of the spin dependent nuclear matrix elements relevant for dark 
matter searches with a variety of experimentally interesting targets (see recent work [8fjl l8lj on such an analysis) . We 
focused on the isovector part of the axial current, which is the most important for the spin induced cross section . 

In particular, the effect on the suitably normalized spin structure functions Ftj of the nucleon form factor has been 
examined. We find that their isospin behavior of the spin structure functions Fij is no longer universal. The isospin 
one modes fall a bit faster as a function of the energy transfer. This effect, however, is not very significant in the 
energy transfer of interest to dark matter searches (see Fig. [4j. In a recent analysis [82|], in which it is shown that it 
is possible to extract all three nucleon cross sections from the data (the coherent as well as the proton and neutron 
spin cross sections), the equality of the spin structure functions is important. In the case of the light nucleus 19 F, 
in which the spin dependent rate maybe more important, we find that even the effect on the integrated structure 
functions Jy, obtained by integrating from zero to the maximum allowed energy transfer, is small: 

Roi = ^-= 0.965948, R n = = 0.936502 
Zoo Zoo 

The effect will be smaller, if the additional u dependence arising from the velocity distribution, not discussed in this 
work, is included. This isospin independent function is also a decreasing function of u, but its precise form depends 
on the WIMP mass. 

In addition to the standard one body contribution, we examined the leading long-range two-body currents, which 
are predicted in chiral EFT [7(J. We confirm the reduction of the isovector matrix elements, quenching factor, by 
as large as 15 — 25%. The precise value is dependent on model parameters, not precisely known. We have found, 
however, that this affects only the static spin ME and it is independent of the target nucleus. Thus this effect can 
be absorbed as a multiplicative factor either in the isovector nuclear matrix elements or better still in the effective 
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coupling parameters, which depend on the assumed particle model, which anyway are going to be derived from the 
experimental data, if and when they become available. In summary the standard analysis of the data previously 
employed is not seriously affected by such effects. 



VI. APPENDIX 



Since the ground state of 19 F is the J = 1/2+ then only the L = 1 multipole contributes (Eq. ([55)) ). In our 
treatment we will separate the couplings from the nuclear operator. As consequence the multi-particle reduced matrix 
elements of the £% =1 and T£ l = \ operators can be written as 

< J\\C\ =1 \\J >= a < J\\£ 5 £=x\\J > +oi < JH^iiiH J > (46) 

where 

< J\\C 6 £i\\J >= < AUT^IIp > + < AIIT^Hp >) (47) 



1 1 1 

< 



J||4iiP >= -^-Wifa) £M J XV2 < X\\T^\\p > + < X\\T^\\p >) (48) 

V V Ap 

with 

/!(<?) = 1 - 6 - - 2 ^ff 2 2c 3 i~ 2T 4—2 ( 49 ) 

We also write 

< J\\T£L\\\J >= a Q < J\\Tll[ m \\J > +a a < J\\Tl l = f 1] \\J > (50) 



with 



and 



with 



< JWULT || J >= ^-^5> A/) (J)(- < A||T( 2 ^)||p > +V2 < AHT^ 1 )^ >) (51) 

X '* X ' Ap 



< J\\TtY\\J >= ±±-j=Mq) ^b Xp (J)(- < AHT^Hp > W2 < AUT^IIp >) (52) 



f 2 ( q ) = l-5-2^ (53) 

M 

The indices A and p run over the single particle orbits of the chosen model space. According to Ref. [29[ this model 
space is the sd one. The quantities a\ p and b\ p are essentially products of the one body coefficients of fractional 
parentage (CFP) for the isoscalar and isovector part of the operator respectively. They depend, of course, on the 
specific interaction and the model space used (see Ref. (29j). 

As it is known the reduced matrix elements of the T^ J ' = \fA~Tiji(qr)\Y l (v) x erp operator is written as 

r-OO 

T (U) =A m R nxl Jr) 3l (qr)R nplp (rydr (54) 
Jo 

where 

r i x 1/2 h ) / , , , \ 

4f = (-l) l H x l p j x j p ljl l p 1/2 j\ Uq Q p J ^6 (55) 
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Therefore the matrix elements (|4"T)) . (|48|) . tfSTT) and (j52")l are written as 

< JII^H J >= ^\^(^'%) + T o (0,1) C 

< J\\^U\J >= ^^AW(V2if J V) + T^\u) 



< J\\T[ l ™\\J >= 



< JHTitf ||J >= -Li * / 2 (u)( - T^\u) + v^lf ^)) 



rf%) + V2T ^( U ) 



732 V4^F 



with 



(0,1)/ 



r, 



( u ) 
o (2,1) W=[( 



1 /2 



6 \5 
1 



Ai + A 2 ) iT - - ( Ax + A 2 ) u + Ai + A 2 



-u/2 



1 2 



15 
ri /2 



"'■ + iVi^)" 2+ (ik A '' 



2 /2 



A' )u 



Tr\u)=[^B 1+ B 2 y~ Z -(B 1 + B 2 ) 



3 V 5 

B 1 + B 2 )u + B 1 +B 2 



-u/2 



1 K ' W 15 1 6 V 5 2 



1 2 



The To' ,J ^(u) part of the operator represents the isoscalar part while Tl l ' J '(u) the isovector one. Using the results of 
our previous shell model calculation for 19 F |29| the coefficients Ax,A 2 ,Bx,B 2 ,A[,A 2 , B[ and B' 2 have been taken the 
values 



1 B [--J-B') U 
15 1 3 V 5 2 J . 

-•{1,3), 



-u/2 



,-u/2 



(56) 
(57) 
(58) 
(59) 

(60) 

(61) 
(62) 

(63) 



A x = 1.19407, A 2 = 1.09082, B x = 1.3695, B 2 = 0.999 



(64) 



and 



A\ = -0.004, A' 2 = -0.0587, B[ = 0.0774, B' 2 = 0.06269 (65) 

The above operators are given in terms of the variable u = q 2 b 2 /2, where q is the momentum transfer and b the 
harmonic-oscillator length. 

According to the above matrix elements the structure functions are written 

Soo = Sn< J||4=il^> 2 + < J\\TlJx°\\J > 2 



(^(«)) a + (r ^(«)) ! 



(66) 



Sxx =8vr < J\\£ 5 L Lx\\J > 2 + < J\\T?Jx L \\J > 



el 5,1 1 



; 2/l 2 ( U ) + / 2 ») (if^cu)) + (a 2 m + 2/Ku)) (if^cu))* 

-! 2V2Tf' 1 )( M )T 1 ^ 1 )( U )(/ 2 ( U )-/ 2 ( M ); 



(67) 



and 



Sox = 8tt< jii^hjx j\\c^x\\J> + <nriLT\\J><J\\V:Jx \\J> 
1 



6 



(2/i(«) + / 2 («))t (2i1) (u)Ti (2,1) («) + V2(/i(«) - /aC*))^ 2 ^)^ ' 1 ^) 
+ V2(/!(u) - Mu))^ ' 1 ^)^ 2 ' 1 ^) + (/ x (u) + 2/ 2 (u))t o (0 ' 1) ( w )T 1 (04) ( u ) 



(68) 
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